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We rapidly change the scattering length a s of a 87 Rb Bose-Einstein condensate by means of a 
Feshbach resonance, simultaneously releasing the condensate from its harmonic trapping potential. 
When a s is changed from positive to negative, the subsequent collapse of the condensate is stabilized 
by the kinetic energy imparted during the release, resulting in a deceleration of the loss rate near 
the resonance. We also observe an increase in the Thomas-Fermi radius, near the resonance, that 
cannot be understood in terms of a simple scaling model. Instead, we describe this behavior using the 
Gross-Pitaevskii equation, including three-body recombination, and hypothesize that the increase 
in cloud radius is due to the formation of concentric shells. 
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Bose-Einstein condensation of neutral atoms is gener- 
ally realized in systems with positive atomic scattering 
length a s ; the resulting repulsive interaction allows the 
formation of large ~ 10 6 atom condensates in harmonic 
potentials. For weakly attractive interactions, however, 
the zero point kinetic energy of the trap can stabilize 
quantum degenerate gases against collapse at sufficiently 
low density [TH3]. Strongly attractive condensates have 
been produced in an important class of experiments that 
uses Feshbach resonances to rapidly tune the scattering 
length a s from positive to negative, but are dramatically 
unstable, resulting in collapsing clouds that expel atoms 
in bursts [1H5]. These works dealt with low density con- 
densates in weak harmonic traps, thereby avoiding strong 
three body recombination that scales as density cubed. 
Here we investigate the stability of a 87 Rb Bose-Einstein 
condensate tuned to negative scattering length in the 
vicinity of a Feshbach resonance. 

In our experiment, we minimize three body recombi- 
nation by simultaneously releasing the atoms from the 
harmonic trap and ramping the scattering length a s to a 
negative value, as done by Volz et al. in Rcf. 8, On the 
high-field side of the resonance, a s is expected to cross 
through at B f — Bq = AB, where Bf is the final mag- 
netic field, Bq is the resonance field, and AB is the reso- 
nance width. For Bf — Bq > AB t h, Volz et al. observed 
the expected decrease in the condensate's Thomas-Fermi 
radius ttf as a function of decreasing Bf, corresponding 
to a decrease in positive a s . For final field settings below 
the expected zero crossing for a s , however, they observed 
an increase in ttf, relative to the minimum cloud width 
observed near the expected zero crossing. They inter- 
preted the swelling in ttf as arising from a destabiliza- 
tion of the condensate in the attractive regime, similar to 
the instability responsible for the explosive condensates 
of Refs. HHS1 

In this work, we observe a similar elevation in ttf of an 
87 Rb condensate, but only for < B } - B < 11(1) f/T, 
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35% lower than the expected crossover to negative scat- 
tering length at Bf — Bq — AB t h — 17 fiT. In contrast to 
the interpretation of Ref. we believe that the observed 
increase in cloud width, rather than indicating instability 
or explosion of the condensate, is related to the sudden 
release of the condensate from the trap, which happens 
much more quickly than the ~1 ms field ramp to neg- 
ative scattering length. There is a short period of time 
after release, thereore, when the mean field energy of the 
condensate is still positive. This positive mean field en- 
ergy is converted to kinetic energy during the brief time 
after release before the scattering length becomes neg- 
ative. In analogy to the the stabilization of a weakly 
attractive condensate by its zero-point kinetic energy in 
a harmonic trap, we suggest that a similar stabilization 
arises from the conversion of mean field to kinetic energy 
prior to the reversal of the mean field from repulsive to 
attractive. In the second case, however, the dynamics 
are more complex, since the condensate, at the moment 
of release, becomes an outward traveling, approximately 
spherical matter wave. Upon reversal of the scattering 
length, it then becomes a complicated superposition of 
inward traveling and outward traveling spherical waves 
that interfere. While we do not directly observe interfer- 
ence phenomena in this experiment, we show that aspects 
of this hypothesis can be successfully modeled using semi- 
analytical arguments as well as numerical simulations. 

We prepare a BEC of N « 2 x 10 6 atoms in a 
crossed optical dipole trap [S], formed by a pair of 
1064 nm laser beams crossing in the x — y plane, with 
final harmonic trapping frequencies of {u) x ,Lu y ,u) z }/2ir = 
{70,55,73} Hz. Our \F = l,m F = -1) BEC starts 
in a small B ~ 1 G bias field along z (vertical) before 
we transfer the atoms to \F — 1, mp — +1) by radio- 
frequency adiabatic rapid passage. (No visible popu- 
lation remains in the \F = l,m F = —1) after Stcrn- 
Gerlach separation of the spin states.) Following the 
transfer, the bias field is set to B « 100.7 mT, near the 
Feshbach resonance under study. 

Because the widest 87 Rb Feshbach resonance has a the- 
oretical width AB t h of only 17 /j,T (170 mG), centered 
at B = 100.7 mT pj |S], tunability of the scattering 
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length requires magnetic field resolution and stability on 
the level of 10 ppm, which we achieve. Our "Feshbach 
coils" are mounted in a Helmholz configuration, and each 
consists of 48 turns of hollow copper tubing. Four sets of 
inlets and outlets provide approximately 5 1/min of cool- 
ing water to each coil. A 300 A current generates the re- 
quired bias field, which settles to the desired set-point in 
less than 1 s [TU] • An additional set of four-turn trim coils 
can rapidly (~ 1 ms) tune the field within ±2 x 10~ 4 T 
(2 G) of the resonance. When current is switched into 
the Feshbach coils, the field overshoots the resonance at 
a high slew rate so that few atoms are lost. Since some 
losses and heating are unavoidable, we allow an addi- 
tional 2 s of free evaporation after the current stabilizes, 
resulting in a nearly pure BEC of « 500 x 10 3 atoms, at 
a field that is either w 100.54 mT or w 100.94 mT. 

For the data of Fig. [T] after the BEC was prepared at 
a field either slightly above or below the Feshbach reso- 
nance, the atom trap was turned off and the trim coils 
were simultaneously ramped by ~ 100 /iT to a final field 
nearer the center of the resonance. The transient field 
from this 500 /is ramp couples inductively to the 48 turn 
Feshbach coils and is sufficient to perturb the power sup- 
ply on a level of ~ 100 ppm. We therefore actively shield 
the larger Feshbach coils commercial supply using a sec- 
ond set of trim coils that have a greater inductive cou- 
pling to the Feshbach coils. The measured transient field 
response associated with a step change of the trim coils 
in the presence of the shielding system has an exponen- 
tial time constant of approximately 3.6 ms. But because 
of the active shielding, the amplitude of the transient is 
much reduced. Following a 0.2 mT (2 G) step change, we 
are therefore able to stabilize the magnetic field to within 
1 /iT (10 mG) on a timescale of ~ 1 ms. This stability is 
maintained for longer than 10 s following the final trim 
coil adjustment. Following release of the trap, and the fi- 
nal 500 /us field ramp, the condensate is allowed to freely 
evolve for 20.0 ms under the influence of the Feshbach 
resonance. The field is then rapidly lowered, almost to 
zero, and the cloud expands for another 9.6 ms. This 
final expansion increases the cloud size and reduces the 
optical depth, facilitating absorption imaging (along z) 
to determine the 2D column density after 30.1 ms total 
time-of- flight (TOF). We fit the 2D image to the sum 
of Thomas-Fermi and Gaussian distributions to obtain 
condensate and thermal characteristics. 

Figure [l] shows the number of Bose-condensed atoms 
Nbec m the condensate along with the Thomas-Fermi 
radius Rtf as a function of the final field setpoint Bf, 
relative to the center of the resonance Bq. Accelerated 3- 
body recombination near the resonance gives rise to the 
sharp loss feature in atom number in Fig. [ija). The loss 
feature is significantly sharper for measurements on the 
low field side (squares) than on the high field side (cir- 
cles). This is expected, since three-body recombination 
scales as p 3 , and much higher cloud densities are expected 
on the high field side where a s is only weakly repulsive 
or even attractive. The loss feature therefore extends to 
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FIG. 1. (color online) (a) The remaining number of condensed 
atoms versus magnetic field reveals a sharp loss feature at the 
resonance field Bo- Data are shown for measurements that 
approach Bo from below (squares) and from above (circles), 
(b) The cloud width after TOF (symbols) is modified by the 
increasing (decreasing) scattering length below (above) the 
resonance. Zero free parameter fits to the data in both (a) 
and (b) are based on both Castin-Dum scaling (solid lines), 
and ID Gross-Pitaevskii simulations (dotted lines). 



higher field offset on the high field side. Somewhat sur- 
prising, however, is the change in curvature of the loss 
feature on the high field side. We naively expect the cur- 
vature of the atom loss curve to remain negative until 
no condensed atoms remain in the cloud. Instead, we 
see a sharp inflection point at Bf — Bq rs 12 /iT, below 
which the slope of the remaining atom number versus fi- 
nal field setpoint decreases significantly. This inflection 
is a general feature observed over a range of experimental 
protocols and suggests a stabilizing influence on the con- 
densate in a region where instability (either rapid collapse 
or explosion) might have been expected. This stabilizing 
influence is aided in part by the losses themselves, which 
decrease the magnitude of the attractive mean field en- 
ergy. 

The Thomas-Fermi radius Rtf [squares in Fig. fTKb)] 
increases dramatically on the low field side of the reso- 
nance by up to a factor of 2 relative to its background size 
of w 100 fim. Although we observed much larger cloud 
widths for fields even closer to the resonance, the clouds 
lose their bimodal appearance, and lack the expected as- 
pect ratio given the initial trap anisotropy [IT]. These 
data points are therefore excluded from Fig. [T] We rely 
primarily on the low field divergence of Rtf to identify 
the center of the resonance, for which radio-frequency 
spectroscopy gives Bq = 100.742(1) mT, consistent with 
previous measurements [5]. 

Approaching the resonance from the high field side, 
the cloud width [filled circles in Fig. [ljb)] decreases as 
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FIG. 2. (color online) The loss rate coefficient K3 can be 
extracted from atom number versus holding time data us- 
ing the analytical model of Eq. |l]), with the prefactor of 
Eq. Q, although this neglects the settling time for a s . Be- 
low the resonance (squares), density is relatively low, and 
the results of this procedure are indistinguishable from a nu- 
merical model based on the Castin-Dum scaling parameters. 
Error bars for the low field data are smaller than the symbols. 
The fit (solid curve) to a power law K3 oc \B; — Bq\^ yields 
Plow = 3.1. Above the resonance, the results of the analytical 
model (open circles) differ from the results of the numerical 
model (closed circles). For data points below 25 /xT, losses 
occur on a timescale of a few ms, faster than the inverse trap 
frequency, and the assumptions of the Castin-Dum numerical 
model break down. A power law fit (solid line) to the high 
field data for Bf — Bo > 25 fiT yields fihigh =3.8. 



a s passes through and becomes increasingly negative. 
We observe a minimum in cloud width at a field offset of 
B - B Q = 11(1) fiT, almost half the value of 20(3) /zT 
obtained in Ref . [5J but comparable to the inflection point 
in atom loss discussed above. We find that the position 
of the field minimum is largely independent of trapping 
frequencies and other experimental parameters such as 
changes to timing protocol. The disagreement between 
our observations and those of Ref. [8J is therefore a mys- 
tery. Comparison to GP simulations (discussed below) 
implies that the position of this minimum occurs several 
pT below the zero crossing of a s , so that it cannot be 
interpreted as a direct measure of resonance width. 

Quantitative analysis of the foregoing atom loss and 
cloud width data requires knowledge of the field depen- 
dence of the 3-body loss rate constant K3 in the vicinity 
of the resonance. Three-body recombination occurs when 
two atoms associate into a molecular state that is deeply 
bound relative to the condensate's ^100 nK tempera- 
ture. Conservation laws require a third atom that gains 
kinetic energy; all three atoms gain sufficient momentum 
to depart the cloud. Three-body losses increase as K 3 p 3 , 
where K 3 is the 3-body loss rate constant and p is the 
local density of atoms within the cloud. As a s decreases 



from positive to negative, the cloud density increases, ac- 
celerating losses. Also, in the vicinity of a Feshbach res- 
onance, the increased overlap with the molecular state 
causes K3 itself to increase dramatically, further increas- 
ing the loss rate. 

K 3 was obtained by measuring the remaining atom 
number N after holding the atoms in the trap at a final 
field Bf for a variable length of time t. These data were 
analyzed in two different ways. First, we fit the resulting 
decay curves to an analytical model 



K X N -K 3 (frn 3 = —K X N - aK 3 N 9/5 , (1) 



dN 

~dt 



where K\ is the one-body loss rate constant and n is the 
density. Two-body losses are assumed to be small [F2] . 
The prefactor a results from the evaluation of the integral 
with the assumption that the BEC retains a Thomas- 
Fermi profile during its evolution, and that a s reaches its 
final value on a time scale that is short compared to the 
losses: 
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Figure [2] shows the field dependence of K3, both be- 
low and above the resonance. On the low-field side of 
the resonance (squares), K3 exhibits a simple power law 
dependence K 3 oc \Bj — -Bo I* 3 with a scaling exponent 
Plow = 3.1(1), slightly higher than the value of approxi- 
mately 2.6 that can be extracted from the data of Ref.fT3l 
In addition to the analytical model of Eqs. |l]) and pj), 
we also fit the data of Fig. [2] to a numerical model based 
on a scaling law [llj , as demonstrated in Ref. [8j with two 
important extensions. First, our numerical solution of 
the Castin-Dum scaling equations incorporates the time 
dependence of a s based on the measured rise-time of the 
field and the predicted field width of the resonance. In 
the vicinity of a Feshbach resonance the scattering length 
scales with magnetic field B as 



a s (t) = a bg 



1 - 



AB 

B{t) - B 



(3) 



where ab g is the background scattering length (w 5.3 nm 
for 87 Rb) and AB and Bq are the width and center of the 
resonance, respectively. We substitutde the theoretical 
value AB t h = 17 p,T [?]. The time dependence of B 
includes a 500 fis ramp to within 1 /iT of the final value, 
followed by an exponential decay to the final setpoint, 
with a time constant of 3.6 ms, as determined by RF 
spectroscopy. 

Our second extension to the analysis of Ref. [8J is the 
self-consistent inclusion of realistic time-dependent losses 
in atom number. We assume a power law dependence 
K3 = a\B — Bq\P , fit the loss data to the Castin-Dum 
scaling equations based on initial guesses for a and j3, 
and obtain a new set of values for K 3 vs B, from which 
new values for a and (3 are obtained. We iterate this 
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procedure until it converges on a stable loss exponent. In 
the case of the low field data, this numerical procedure 
yields results that are equivalent to those obtained from 
the analytical model. 

For the high field data, however, there is a striking dif- 
ference between the results of the analytical model (open 
circles) and the numerical model (filled circles). This 
is because the time dependence of a s , neglected in the 
analytical model, is more important on the high field 
side, where decreasing a s shrinks the Thomas-Fermi ra- 
dius of the cloud, increasing density, and accelerating 
losses. For values of \Bf — Bq\ > 25 /iT, the numer- 
ical model yields results for K 3 that appear to obey a 
power law with Phigh — 3-8(1), significantly higher than 
Plow We note that, based on Eq. ([3]), K 3 also obeys a 
power law with \a/ab g — 1|, which mirrors the abscissa 
of Fig. |2j rather than with \a/a,b g \ directly. The uni- 
versal scaling law K 3 cx a 4 [Uj, which has been con- 
firmed for a I dbg 3> 1 |15| . is therefore not obeyed for our 
data. However, this universality law is not expected to 
to hold for small a, or even very close to the resonance 
[16) , and we are therefore not necessarily surprised to find 

Phigh 7^ Pi 

For the high side data with \B f - B \ < 25 /iT, K 3 
appears to diverge rapidly. This is most likely due to 
a departure of the cloud density from an ideal Thomas- 
Fermi profile, as the loss rate approaches the ~ 1/ui equi- 
librium timescale for the cloud, where uj = (cl^u^u^) 1 / 3 
is the geometric average trapping frequency. Fits to the 
loss data in this regime are visibly poor and yield large 
uncertainties. For subsequent analysis of the TOF data 
of Fig. [TJ we have assumed K 3 cx (B — i? ) 3 ' 8 , even for 
B - B < 25 fiT. 

We now attempt to fit the untrapped behavior of Fig. [I] 
using the same Castin-Dum scaling argument that was 
just used to analyze the trapped loss data of Fig.[2j A zero 
free parameter theory curve (solid line) , using the values 
for Pi ow and Phigh obtained above, reproduces the loss 
feature of Fig. [lja) quite well, including the inflection in 
the slope of the loss curve on the high field side. Based on 
this model, we understand this reduction in loss rate in 
terms of a reduction of the mean field energy that drives 
the collapse, since mean field energy depends on atom 
number. 

The Castin-Dum theory curve of Fig. |T|(b) likewise 
yields reasonable agreement for Bf < Bq, where the max- 
imum cloud width (excluding points closer to Bq that ap- 
pear to have a thermal profile) corresponds to a factor of 
35 increase in a s . On the high field side of the resonance, 
the Castin-Dum model once again produces an inflection 
in the slope of Rtf vs Bf — Bq. As noted previously, 
rapid losses reduce the negative mean field energy that 
drives the collapse, thereby slowing the collapse. How- 
ever, the Castin-Dum model has Rtf without an 
upturn at very low fields, in disagreement with the ex- 
perimentally observed minimum and subsequent increase 
in cloud width with decreasing field offset. Rather, to re- 
produce the observed increase in cloud width, which has 
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FIG. 3. (color online) a) Cloud density is mapped to a color 
scale, showing the radial density profile as a function of time 
during TOF. After an initial expansion, the cloud contracts 
into concentric rings until t — 20 ms when the field is turned 
off. b) A line cut through the image of a) at t = 30 ms shows 
a second large peak in the radial density profile (dotted line) , 
which should be visible in the integrated OD (solid black line) . 
A small modulation is visible in the experimentally observed 
OD profile (open circles), relative to a Thomas-Fermi plus 
Gaussian fit (solid red line). 



previously been interpreted as an instability [5] , we must 
model the system within the Gross-Pitaevskii equation 
(GPE). 

To better understand the observed inflection in loss 
versus field in the negative a s regime, we modeled our 
system within the GPE, which has the form of the nonlin- 
ear Schrodinger equation. Assuming spherical symmetry, 
we cast the 3D GPE as an effective ID radial equation 

mm, 



dt 
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+ V(r,t) 



K 3 (t)\i;(r,t)\ 4 



(4) 



where the time-dependent radial wavefunction ip(r,t) is 
normalized to the total atom number N. The poten- 
tial V(r,t < 0) = mu 2 r 2 /2 while V{r,t > 0) = 0, cor- 
responding to the release of the atoms from the trap 
at t = 0. The time-dependent interaction strength 
g(t) = Airh 2 a s (t)/m where a s (t < 0) = ab g and a s (t > 0) 
is given by Eq. d3| with the inclusion of a time-dependent 
magnetic field B(t). The phenomenological loss term is 
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equivalent to Eq. ([I]) but is completely general, making 
no assumptions about the form of the density profile. 
The loss constants Ki and K 3 (B(t)) are taken from fits 
to the data of Fig. [2] The time-dependence of the field 
is the same as for the Castin-Dum analysis above. 

We numerically solved the GPE using the Crank- 
Nicolson (CN) method jTSJITS] ■ Beginning with a trapped 
condensate (u/2ir = 65 Hz) at a magnetic field of 
B — Bq = 0.2 mT, the ground state was determined by in- 
tializing the wave-function with a Thomas-Fermi profile, 
then propagating in imaginary time to eliminate higher 
order spatial modes [TB], resulting in the true ground 
state. The wavefunction was normalized to N particles 
after each imaginary time-step, where N ranged from 
2.5 x 10 4 to 2.5 x 10 5 atoms. The CN algorithm was 
then run in real time with the trapping potential turned 
off; the field was ramped to a final field Bf, and the 
atom number was allowed to diminish according to the 
loss constants. After t=20 ms, the magnetic field was in- 
stantaneously set to 0. The simulation was then allowed 
to evolve for an additional 10 ms, corresponding to a total 
of 30 ms TOF, in order to provide direct comparison with 
experimental results. Finally, the radial density profiles 
were integrated along one dimension so that the resulting 
profiles correspond to the optical depth (OD) profile that 
is obtained experimentally from an absorption image. 

Fig. [3ja) shows the simulated ID radial density profile 
evolving during time-of- flight (TOF). Following release at 
t = 0, the BEC containing 2.2 x 10 5 atoms undergoes ini- 
tial expansion for ~ 1 ms before the field settles to a value 
of B — Bq = 10 /iT, where a s — — 0.7a& g . Despite the 
negative scattering length, corresponding to attractive 
mean field energy, the condensate does not immediately 
collapse. Rather, the kinetic energy imparted to the sys- 
tem during the brief period following release (while the 
mean field was still positive) imposes a quadratic spatial 
dependence on the phase of the wavefunction. Subse- 
quently, when the wavefunction begins to be pulled in- 
ward by attractive mean field energy, the non-uniform 
spatial dependence of the BEC density reverses the mo- 
tion of outward going motion in a non-uniform way, caus- 
ing the wavefunction to interfere with itself, rather than 
to uniformly collapse. This explains the arms that begin 
to form in the ID radial density map around t = 5 ms. 
Translated into 3 dimensions, these arms correspond to 
the development of concentric shells. At t = 20 ms, the 
field is turned off completely, and the scattering length 
becomes positive once again. The concentric shells ex- 
pand and blur, but are still visible in the simulated image 
at 30ms TOF. 

Fig. [3jb) shows a line cut through the radial density 
image ofFig.|3];a)atTOF = 30 ms. The simulated ID ra- 
dial density profile (dotted line) shows a strong secondary 
peak corresponding to a shell of atoms at r « 20 /im. 
The optical depth profile (solid black line) is calculated 
from the 3D density distribution corresponding to the 
radial density profile, by integrating along a Cartesian 
axis. The prominent secondary peak in the radial den- 



sity profile is still visible in the calculated OD profile, 
which is what should be observed in experiment. A line 
cut from the experimentally observed OD profile (open 
circles) shows a small amount of spatial modulation com- 
pared to the Thomas Fermi plus Gaussian fit to the data 
(solid red line), but similar artifacts can sometimes be 
seen far from the Feshbach resonance, and are decidedly 
inconclusive. Note that the a;-scale for the experimental 
data differs from the r-scale for the simulation results. 

The CP simulations were performed at several differ- 
ent fields to obtain final cloud characteristics for com- 
parison to Fig. [T] For small B — Bq, the appearance of 
the secondary peak presents a challenge to identification 
of the Thomas-Fermi radius. Our procedure is to ignore 
the central peak which can, in simulation, become quite 
prominent, despite the inclusion of realistic 3-body losses. 
We mask the data to select only the monotonic region of 
the secondary peak and fit this to a Thomas-Fermi (in- 
verted parabola) profile. Fig. [I] shows the results of the 
Gross-Pitaevskii (GP) simulation (dotted lines) for the 
high side of the resonance. As for the CD numerical solu- 
tion, the GP simulation reproduces the number loss data 
quite well, including the inflection in the slope of the loss 
curve. In Fig. [ijb), the GP simulation underestimates 
the magnitude of Rtf, but qualitatively reproduces the 
upturn in Rtf at low values of Bf — B . 

The upturn in Rtf is not observed in our experiment 
at Bf — Bq = AB ps 17 /iT, as previously reported in 
Ref. HI but rather at a much lower field, less than 10 /iT. 
This difference suggests that the upturn in Rtf is a poor 
measure of resonance width AB which, to our knowledge, 
has otherwise never been measured. We therefore con- 
sider whether there are any other prominent features in 
the observed resonance data, particularly near the ex- 
pected crossover from positive to negative a s , which oc- 
curs at B f — Bo = AB. As shown in Fig. |2Fa), we have 
measured atom number versus offset field on the high 
side of the resonance, for several clouds with varying ini- 
tial atom number. This is accomplished by preparing 
the cloud as for the data of Fig. [I] but varying the MOT 
loading time from a few hundred ms to several seconds. 
For the three datasets with highest initial atom number 
Ni = 2.3 x 10 5 , 1.7 x 10 5 , and 1.5 x 10 5 , the observed loss 
feature is unchanged at the high side data of Fig. [TJa) . 
GP simulations (solid lines) show reasonable qualitative 
agreement with these data. For lower N{ — 5.1 x 10 4 and 
3.3 x 10 , the atom number is nearly independent of final 
field offset until Bf — Bq ps 14 /iT, which appears as a 
threshold field for atom loss. The GP simulations, which 
assume AB = 17 /iT, do not reproduce this threshold 
behavior. 

The threshold behavior observed in Nbec is accom- 
panied by an increase in thermal atom number N t h erm ai j 
as shown in Fig. [4jb) . The increase in thermal number 
is abrupt for all values of Ni, and occurs at a final offset 
field Bf — Bo fa 12 /iT that is significantly lower than 
the theoretical width AB = 17 /zT, which is marked by a 
dotted line in the figure. The presence of thermal atoms 
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FIG. 4. (color online) a) Number of condensed atoms Nbec 
measured after TOF versus final offset field Bf ~ Bo- For 
low initial atom number Ni — 51 k and 33 k atoms, a loss 
threshold is visible at B — Bo ~ 14 /iT, significantly less that 
the theoretical resonance width AB — 17 /iT, which is in- 
dicated by the vertical dotted line. Solid lines indicate GP 
simulation results, which assume AB = 17 /iT and which 
fail to reproduce the observed loss threshold for Ni — 51 k 
and 33 k atoms, b) Number of thermal atoms Nrhermai after 
TOF. A threshold for heating is visible at B f - B « 12 fj,T. 
c) Thomas-Fermi radius Rtf vs final offset field. The mini- 
mum in J?tf occurs well below the expected zero crossover for 
a s , both for experimental data (symbols) and GP simulations 
(solid lines). 



trap symmetry, the blurring that occurs in the final 10 ms 
of expansion, a lack of spatial resolution in our imaging 
setup, or "seeding" of the gain processes by thermally 
induced modulations in density. However, the formation 
of this shell structure can be viewed as an interference 
effect, and self- interference has been associated with tur- 
bulence, vortex formation, and heating |20| . It is con- 
ceivable that the onset of heating in our experiment is an 
indication of self-interference of the condensate as it col- 
lapses under the influence of attractive interactions. In 
simulations, the onset of interference effects, manifest in 
the formation of a secondary peak in the density profile, 
occurs within 1 — 2 /iT of the zero crossing of a s . We 
therefore have a second possible indication of the true 
resonance width. 

Finally, we see in Fig. [3Jc) that the increase in Rtf 
at very low values of Bf — Bq < 10 /iT is a robust phe- 
nomenon for any Ni. We interpret this phenomenon in 
terms of the shell structure that appears in the density 
distribution of Fig. [3] Interference modifies the usual 
Thomas-Fermi (inverted parabola) density distribution, 
such that the contracting cloud of atoms leaves behind 
a shell of atoms at a radius that is larger than a sim- 
ple Castin-Dum scaling law would predict. Here we have 
to posit a conspiracy of heating and other imperfections 
that transform this shell of atoms into an effective in- 
crease in the overall diameter of the cloud. This is clearly 
unsatisfying. We note that the GP simulations in this pa- 
per solve the 3D GPE as an effective ID radial equation. 
It is possible that a full 3D simulation including initial 
trap assymetries would fail to produce the shell structure 
of Fig. [3] or would transform it in some other way. 

We have shown that the collapse of an untrapped 
87 Rb condensate with negative scattering length can be 
made to proceed in a stable manner. The condensate 
is stabilized against explosive collapse in part by the ki- 
netic energy imparted upon its release from the trap. 
The observed increase in Rtf may be related to a self- 
interference effect that imposes a shell structure onto 
the condensate, as seen in Gross-Pit aevskii simulations. 
While this shell structure is not observed directly in ex- 
periment, self-interference may be responsible for turbu- 
lence resulting in the observed increase in thermal frac- 
tion at low offset field. 



is not included in the GP simulations, since the 3-body 
loss mechanism is usually associated with loss of all three 
atoms from the trap. On the low side of the resonance, 
the existence of a weakly bound state could mediate three 
body recombination at sufficiently low energies that the 
atoms remain trapped. On the high side of the resonance, 
however, no such weakly bound state exists. 

We recall that the shell structure of Fig. [3] is not ob- 
served in our experiment. This may be due to imperfect 
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